New Universality Class at the Magnetic Field Tuned Superconductor-Insulator 

Transition? 
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The superconductor-insulator transition in ultrathin films 
of amorphous Bi was tuned by an applied magnetic field. A 
finite size scaling analysis has been carried out to determine 
the critical exponent product vz — 0.7 ±0.1 for five films of 
different thicknesses and normal state resistances. This result 
differs from the exponents found in previous experiments on 
InO x and MoGe films. We discuss several possible reasons for 
this disagreement. 

PACS numbers: 74.76.-w, 74.40. +k, 72.15.Rn 



The superconductor-insulator (SI) transition in two- 
dimensional disordered systems has been extensively 
studied theoretically in the context of the dirty boson 
model at zero temperature. Monte Carlo simulations 
jj],D and renormalization group calculations || mostly fo- 
cused on the case when the transition is tuned by chang- 
ing the carrier concentration, coupling strength or the 
strength of the disorder. The SI transition can also be 
tuned by changing the external magnetic field Scal- 
ing arguments [^) give a lower bound on the correlation 
length exponent v > 1, and predict the value of the dy- 
namical critical exponent to be z=l An experimen- 
tal scaling analysis of the magnetic field tuned SI tran- 
sition has been carried out for thin films of InO x JtJ and 
MoGe ||. Both experiments found vz w 1.3, consistent 
with the theoretical predictions. We have studied the 
magnetic field tuned SI transition in amorphous ultra- 
thin homogeneous films of Bi. Finite size scaling analysis 
yields vz « 0.7 which is in disagreement with the predic- 
tions of Ref. [Q, as well as with the results of Refs. |7]] 
and§. 

The investigations were carried out on five ultrathin 
Bi films with thicknesses close to 12-13A. The films 
were evaporated on top of a 10A thick layer of amor- 
phous Ge, which was pre-deposited onto a 0.75mm thick 
single-crystal of SrTiO% (100). The film thickness was 
systematically increased in situ under UHV conditions 
(~ 10~ 10 Torr), with the substrate temperature kept be- 
low 20 K during all depositions. Film thicknesses were 
determined using a previously calibrated quartz crystal 
monitor. Films prepared in this manner are believed to 
be homogeneous, since it has been found that they be- 
come connected at an average thickness on the order of 
one monolayer Between the successive depositions, 
resistance measurements were carried out using a stan- 



dard dc four-probe technique with current bias up to 
50nA. Magnetic fields up to 12kG perpendicular to the 
plane of the sample were applied using a superconducting 
split-coil magnet. 

A typical temperature dependence of the resistance as 
the magnetic field changes is shown on Fig.|l|. In zero 
field, the resistance decreases with increasing tempera- 
ture suggesting the existence of superconducting fluctua- 
tions. Magnetic field destroys this downward curvature, 
and at some critical magnetic field, B c , the resistance is 
independent of temperature. In magnetic fields higher 
than B c the film is insulating, with 5R/5T > 0. The 
critical magnetic field is determined by plotting the re- 
sistance as a function of magnetic field for different tem- 
peratures (Inset of Fig.^) and identifying the crossing 
point for which the resistance is temperature indepen- 
dent, or by plotting dR/dT as a function of magnetic 
field at the lowest temperatures and finding the field for 
which (dR/dT) =0@. 

The scaling relation for the resistance of a two dimen- 
sional system close to the magnetic field tuned supercon- 
ductor insulator transition was suggested by Fisher M: 



R(5,T) = R C /(ST- 1 ^*) 



(1) 



Here 5 = B — B c is the deviation from the critical field, 
R c is the critical resistance at B = B c , f(x) is a universal 
scaling function such that /(0) = 1, v is the coherence 
length exponent, and z is the dynamical critical expo- 
nent. To carry out the analysis of the data, we introduce 
a parameter t = T~ x l vz , and adjust t(T) at each temper- 
ature so that all the resistance data collapse on a single 
curve. The critical exponent product vz is then found 
from the temperature dependence of t. This procedure 
does not require prior knowledge of the critical exponents 
or detailed knowledge of the functional form of the tem- 
perature or field dependence of the resistance. 

The collapse of the resistance data as a function of St 
for one of the samples is shown in Fig.||. The critical 
exponent product vz, determined from the temperature 
dependence of the parameter t (Inset of Fig.||) , is found to 
be vz = 0.7 ±0.1. The same exponents were obtained us- 
ing a somewhat different method of plotting (dR/ dB)\s c 
vs T _1 on a log-log plot and determining the slope which 
is equal to l/vz [Q. The analysis was carried out for five 
films of different thicknesses, all yielding vz = 0.7 ± 0.1, 
apparently independent of the film thickness. In fact, 
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plotting the normalized resistance R/R c as a function 
of the scaling variable reveals that they all scale together 
(Fig. ||). This behavior strongly suggests that the scaling 
function is universal. 

The critical resistance, however, does not seem to be 
universal. Figure || shows that R c decreases as the critical 
field increases, roughly in a linear fashion. Very similar 
behavior was observed by Yazdani and Kapitulnik J^] . 

The critical exponent product vz found in this experi- 
ment does not agree with the theoretical prediction that 
v > 1 and z = 1 Q for a disordered system, nor with the 
experimentally obtained z — 1 and v w 1.3 for a mag- 
netic field tuned SI transition in thin films of amorphous 
InO x @ and MoGe J§. There could be a number of 
possible reasons for this disagreement, which we discuss 
next. 

The prediction that v > 1 follows from an exact the- 
orem |To| which is expected to hold for any transition 
which can be tuned by changing the strength of the dis- 
order. It has been argued recently that this is not 
necessarily the case, and that the disorder averaging may 
play an important role in determining the critical expo- 
nents. The nature of disorder would be crucial here. 

Without disorder, the SI transition in a 2D system 
might belong to the universality class of the classical 3D 
XY model The exponent product we find for the 

field-driven transition is very close to that obtained for 
the classical 3D XY model. Numerical simulations of a 
T=0 critical point for the boson Hubbard model in 2D for 
which there is no disorder Q also find vz ss 0.7. Since our 
films are believed to be very homogeneous, it is possible 
that the length scale which characterizes the disorder is 
very small compared to the diverging coherence length, 
and the system behaves as if there were no disorder. 

The importance of the length scale and the nature of 
disorder might also explain the difference between our re- 
sults and those of Refs. Jt| and Since both InO x and 
MoGe are composite materials, the length scale for dis- 
order might be quite different from that of monoatomic 
amorphous Bi. Indeed, InO x films in a magnetic field 
appear to be in a quasi-reentrant regime, suggestive of 
the behavior of granular films in zero field |12j| . Second, 
the way the samples are prepared might result in a pro- 
foundly different nature of the disorder. In InO x and 
MoGe, the normal state resistance is determined by the 
composition (content of oxygen and germanium, respec- 
tively), while in quench condensed Bi films it depends on 
the film thickness. We also wish to point out that the Bi 
films used in this study are almost an order of magnitude 
thinner than the samples of Refs. j?]] and ||. 

It has been suggested that a local dissipation for the 
phase of the superconducting order parameter due to 
gapless electronic excitations might change the universal- 
ity class of the system and lead to a non-universal critical 
resistance |13 . The critical resistance is expected to in- 



crease with increasing damping due to dissipation, which 



might be expected to increase with decreasing normal 
state resistance. We however observe that the critical 
resistance decreases as the normal state resistance goes 
down, which is exactly the opposite of the behavior pre- 
dicted by Ref. @. 

A recent work by Shimshoni et al. |l4j considers the 
system close to the SI transition to be a binary compos- 
ite of an insulating and a superconducting phase. Using 
incoherent Boltzman transport theory, they derive resis- 
tivity laws in different temperature regimes and predict 
finite dissipation at T=0 at all fields. We do not observe 
any saturation in the temperature dependence of the re- 
sistance as the temperature decreases (8R/ST is non-zero 
down to the lowest temperatures). A satisfactory fit to 
the predicted resistivity laws also could not be obtained. 

Local fluctuations of the amplitude of the supercon- 
ducting order parameter are usually neglected in the 
scaling theory and the numerical simulations, but might 
actually be important. It has been suggested that the 
percolation of islands with strong amplitude fluctuations 
might mask the localization exponents obtained from the 
scaling theory JTsf . 

The existence of electronic excitations was also used 
to explain the non-universal behavior of the critical re- 
sistance in Ref. (8). A model of two-channel conduction 
was proposed, where the conductance due to the elec- 
tron channel adds to the conductance due to the boson 
channel. Films with lower normal state resistances would 
have lower critical resistances due to a larger fraction of 
normal electrons. Our data could only be explained this 
way if the quantum resistance due to pairs were greater 
than h/Ae 2 . The conductance due to the electronic chan- 
nel in a magnetic field might also depend on the strength 
of the spin-orbit interactions, which is another major dif- 
ference between our samples and those of Refs. M and 
p|. There is, however, a striking similarity in the mag- 
netic field and normal state resistance behavior of the 
critical resistance. 

We should note that the critical resistance is predicted 
to be universal only at T=0, while the finite temperature 
corrections are expected to be scaled by T c B . Normal 
state resistances of the MoGe films of Ref. M are a fac- 
tor of 3-10 lower than the Bi films considered here, and 
are far away from the critical point of the disorder driven 
transition. This means that our samples are probing a 
different part of the phase diagram @Q, and the finite 
temperature corrections might be more important in one 
case then the other. There is also a question of whether 
any of the experiments really accessed the quantum crit- 
ical region |l^| . 

In order to resolve some of the issues discussed in this 
paper, simulations in which the magnetic field is the tun- 
ing parameter would be very useful. 

We gratefully acknowledge useful discussions with A. 
P. Young, S. Sachdev and S. L. Sondhi. This work was 
supported in part by the National Science Foundation 
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FIG. 1. Resistance per square as a function of tempera- 
ture in different magnetic fields, ranging from OkG (bottom) 
to 12kG (top), with IkG increments. Inset: Resistance per 
square as a function of magnetic field at 0.15, 0.17, 0.19, 0.2, 
0.25, 0.3 and 0.35K. 
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FIG. 2. Resistance per square as a function of the scal- 
ing variable, t\B — B c \, for seventeen different temperatures, 
ranging from 0.14K to 0.5K. Here t = T~ 1/vz is treated as 
an adjustable parameter to obtain the best collapse of the 
data. Different symbols represent different temperatures. In- 
set: Fitting the temperature dependence of the parameter t 
to a power law determines the value of vz. 

FIG. 3. Normalized resistance per square as a function of 
the scaling variable, T- 1/l,z \B - B c \, at 0.15K for five films: 
12.239A (squares), 12.353A (diamonds), 12.455A (crosses), 
12.553A (triangles) and 12.660A (circles). Only one temper- 
ature was shown for each film for clarity. 

FIG. 4. The critical resistance as a function of the critical 
field for the same five films of Fig. 3. Here R c decreases with 
increasing thickness, as thicker films have lower normal state 
resistances and higher critical fields. 
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